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Abstract. The absolute integrated intensity of the singlc-phonon Raman peak 
at 1580cm -1 is calculated for a clean graphene monolayer. The resulting intensity 
is determined by the trigonal warping of the electronic bands and the anisotropy 
of the electron-phonon coupling, and is proportional to the second power of 
the excitation frequency. The main contribution to the process comes from 
the intermediate electron-hole states with typical energies of the order of the 
excitation frequency, contrary to what has been reported earlier. This occurs 
because of strong cancellations between different terms of the perturbation theory, 
analogous to Ward identities in quantum electrodynamics. 
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1. Introduction 

In the past decades, Raman spectroscopy [I] techniques were successfully applied to 
carbon compounds, such as graphite [2] and carbon nanotubes [3]. Upon the discovery 
of graphene [1] , Raman spectroscopy has proven to be a powerful and non-destructive 
tool to identify the number of layers, doping, disorder, strain, and to characterize the 
phonons and electron-phonon coupling [5] [T] [9] [101 E] • 

The most robust feature in the Raman spectra of all conjugate carbon compounds 
is the so-called G peak, corresponding to the in-plane bond-stretching optical vibration 
of sp 2 -hybridized carbon atoms. In graphene this vibration has the E 2g symmetry, it is 
doubly degenerate, and its frequency iv p h ~ 1580 cm -1 . While the G peak frequency 
is sensitive to external factors, such as doping level [6] [T] [8] or strain [121 EH HI], 
its total frequency-integrated intensity Iq is often assumed to remain constant under 
the change of many external parameters, depending only on the excitation frequency. 
Thanks to this robustness, Iq is often used as a reference to which intensities of 
other peaks are compared [T3] Q31 [T3] , since measurement of absolute peak intensities 
represents a hard experimental task (the first use of Iq as a reference for other Raman 
peak intensities in graphite probably dates back to 1970 (16)). 

Given this popular role of Iq as a reference, it is clear that having at hand a 
theoretical expression of the absolute intensity in terms of the basic parameters of the 
material (such as the electronic dispersion, the electron-phonon coupling, etc.), would 
be useful. Nevertheless, to the best of the author's knowledge, no such expression is 
available at present. Iq was discussed in a recent paper by the author [17j . where 
the conclusion was reached that even in the limit when the excitation frequency Lo- m is 
small compared to the energy scale to characterizing the electronic dispersion (several 
eV), the intensity Iq is contributed by the whole conduction and valence bands, not 
just by low-energy states in the vicinities of the Dirac points. As will be shown below, 
this conclusion is wrong, since Ref. [17] missed strong cancellations occurring as a 
consequence of a Ward identity. This affects the dependence of Iq on u>- ln . 

In the present work Iq is calculated for a clean graphene monolayer suspended in 
vacuum. When aj- m <C to, the intensity is indeed determined only by electronic states 
with energies ~ u)- m , not the whole conduction and valence bands (however, these 
energies do not have to be close to wi n /2). In this limit Iq can be expressed in terms 
of a few parameters, characterizing these low-energy states. As correctly mentioned in 
Ref. |17j , trigonal warping of the electronic bands and the anisotropy of the electron- 
hole coupling turn out to be crucial for the G peak. In the limit w p h <C Wj n <C to the 
frequency dependence is Iq oc uif n . At higher frequencies, the calculation is performed 
using the nearest-neighbor tight-binding model. The main feature found is a strong 
enhancement of Iq when the frequency matches the energy of electron-hole separation 
at the M point of the electronic first Brillouin zone, where a van Hove singularity in 
the electronic density of states occurs. 

2. Calculation 

The Bloch form of the electronic wave function in the tight-binding model is 




(1) 
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Figure 1. Honeycomb lattice with the two sublattices A and B and the three 
nearest-neighbor bond vectors di,d2,d3. 



where ipn(r) is the atomic wave function localized near the origin, a = A, B labels 
the sublattices, and the two-dimensional integer vector n labels the unit cells, each 
containing two atoms from the two sublattices. The wave vector k runs over the first 
Brillouin zone. The Bloch amplitudes U a k satisfy the Shrodingcr equation (we neglect 
the non-orthogonality of the atomic orbitals) 



e-iT, 



(0) 



= 0, 



(2) 



where the Hamiltonian matrix is given by 
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Here in ~ 3 eV is the nearest-neighbor coupling matrix element, and di : 2,3 are the 
vectors connecting an A atom to its three nearest neighbors, as shown on Fig. [1] 
Their length is |di,2,3 1 — a ~ 1.42 A. The energy eigenvalues are e = ±ek, where 
ek = |.Hk|- Note that is not a periodic function in the first Brillouin zone; its 
values on the boundaries related by a reciprocal lattice vector differ by a constant 
phase factor e ±27n / 3 . 

Suppose that all atoms belonging to the sublattice a are displaced by u Q , and 
ua = —ub- Such a uniform displacement corresponds to optical phonons with the 
wave vector q = 0. These are the phonons responsible for the G peak, as fixed by the 
momentum conservation (the wave vectors of both incident and scattered photons are 
negligibly small). In the tight-binding model the mechanism of the electron-phonon 
coupling is the change of the nearest-neighbor electronic matrix element to due to the 
change in the bond length a, dto/ 'da ~ — 6eV/A. Then, to linear order in the atomic 
displacements, the electron-phonon interaction Hamiltonian is given by 

u A -u B ( F k \ F 2i dtp dH k 
y ' k tpa da 9k 



V, 



(o) 



(36) 



Calculation of the Raman G peak intensity in monolayer graphene 



4 



The electromagnetic field of the incident and scattered light is described by the 
vector potential A. The Hamiltonian describing interaction of the electrons with the 
field can be obtained by the Peierls substitution k — > k — (e/c)A in the electronic 
Hamiltonian. Since we are interested in a two-photon process, we should expand the 
Hamiltonian to the second order in A: 

(3c) 

(3d) 

(3e) 
(3J) 
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Vk ~ dk ' Wk ~ dhdkj ' 5k " 9% ' k " fl^flfe, ■ 1 9) 

Here i = x,y label the in-plane Cartesian components, and the summation over 
repeated indices is assumed. Note that for the description of resonant multiphoton 
processes it is sufficient to keep term only, as it was done in [T7], since the 

contribution of other terms is (i) off-resonant and (ii) smaller by a factor ~ oj- ln /to. 
However, for the one-phonon process studied here the contributions from all terms 
turn out to be of the same order; in fact, they almost cancel each other. 

To calculate the Raman scattering probability, we follow the general steps of 
Ref. [17] : pass to the second quantization, calculate the scattering matrix element, 
and substitute it into the Fermi Golden Rule. The zero-temperature electronic Green's 
function is a 2 x 2 matrix, 

G k (e) = (e-JJk + iO+sgne)- 1 . (4) 

Upon quantization of the phonon field the atomic displacements corresponding to 
the optical phonons are expressed in terms of the phonon creation and annihilation 
operators bL^b^ as 

.00 



UA.n = — U.B • 



V 



(i q ,e iqr °+Lc.) . (5) 



Here fj, = x,y labels the two degenerate phonon modes, e^) is the unit vector in the 
corresponding direction, M is the mass of the carbon atom, N is the total number of 
the carbon atoms in the crystal, and "h.c." stands for the hermitian conjugate. A small 
finite wave vector q has to be introduced in order to treat the in-plane momentum 
conservation properly; only the q = modes will contribute to the final result. The 
transverse electromagnetic field is quantized inside the large volume L x L y L z in the 
Coulomb gauge: 

e (Q/) a Q/ e iQr + h.c.l . (6) 




L x L y L z Q 



Here Q is the three-dimensional wave vector, I — 1,2 labels the two transverse 
polarizations along the two unit vectors e^^^ _L Q. 
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The probability for an incident photon with the polarization e ln and frequency 
Wj n to be scattered into an element of the solid angle do out with the polarization e out , 
and with the frequency uj out = u>- m — oj pn fixed by the energy conservation, is given by 

^ G - ^out |o KAijl pine „out\* (/*) 2 t 7 \ 

d 0out - (2iry c *^r M 1 [6j > e i ■ [n 

Here M. % ^ 1 is the transition matrix element (the factor of 2 explicitly takes care of 
the two spin projections of the electron). It is given by the sum of the following 
contributions: 



M ijl = ^e 2 / L x L y r d 2 k ale ^ 



^/^in^out y 2NMtv ph J (2ir) 2 2-k 
x Tr{V\ jl + vf + vf + vf + vf + vf } , (8) 

Vf = G k (eRG k (e - w in )F k G k (e - uw)^ (9a) 

vf = G k (e)v J k G k (e + w out )^G k (e + u in )vl (96) 

Vf = G k (e).^G k (e - w out )^, (9c) 

Vf = G k (e)v k G k (e-u >in )g l k j , (9d) 

W =G(e)^G k (e-w ph )F k , (9e) 

Vf = ^Gk(e). (9f) 

H prescribes closing of the e-integration contour in the upper half- 
plane. In fact, it is important only for the last term, Eq. which corresponds to 
the sum over the filled valence band. Upon integration over e, Eq. ([8]) reduces to 
the standard perturbation theory expression for the Raman amplitude as the sum 
over intermediate states. Each intermediate state contains an electron with the wave 
vector k in the conduction band and a hole in the valence band with the wave vector 
— k. The k integration is performed over the first Brillouin zone. The six terms given 
by Eqs. (|9a|) -([9^) can be represented pictorially by the diagrams shown in Fig. [2] 
In fact, V\ k = V\ k because the Green's function and all the vertex matrices satisfy 
G£(e) — a x G k (e)a x and G k (— e) = — a z G k (e)a z (the T superscript denotes the matrix 
transpose, and a x ^ z are the Pauli matrices). 

If one integrates each of the terms in Eqs. (|9al) -([9^) separately, the main 
contribution to the k-integral comes from the "bulk" of the first Brillouin zone, rather 
than the vicinities of the Dirac points (as it was pointed out in Ref. 17 \). Let us 
formally consider, however, the whole expression for M*i at ojj n = oj out = 0. Then 
the sum of the terms (p?a[) -(|9"^) is a total derivative, d 2 [F k G k (e)]/(dkidkj). Thus, 
the k-integral vanishes, i. e., at u>i n — uj ou t — all terms contributing to the matrix 
clement, cancel. This cancellation is far more general than the tight-binding model, 
used here. Indeed, consider linear response of the electronic current j to a static 
homogeneous vector potential A. The corresponding response function Xij must 
vanish, since an observable (current) cannot respond to a pure gaug^J This holds 
for any configuration of the atomic positions, hence its derivative with respect to 

X Strictly speaking, it is the limit lim q _>o lim^^o Xij (q> w ) which vanishes, while the calculation 
done in this work corresponds to the opposite order of limits. The two limits commute for undoped 
graphene, when the Fermi surface has zero length. At finite doping, the contribution of the filled 
valence band states far below the Fermi energy still cancels out. 
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Figure 2. (Color on line.) Diagrams corresponding to Eqs. I|9all -||9/||. Solid 
arrows correspond to the electronic Green's functions, the dashed arrow represents 
the phonon emission vertex, the zigzag arrows represent the photon absorption 
and emission vertices. 

the displacements, dxij/dui, must also vanish. As follows from the Kubo formula, the 
matrix element Ad 1 ^ 1 at w; n = 07 out = is equal just to this derivative, up to a constant 
factor. This fact is analogous to the Ward identity in the quantum electrodynamics. 
The cancellation of different terms contributing to JA lJ was overlooked in Ref. |17j . 
and it was wrongly concluded that M. % ^ 1 is determined by the whole of the first 
Brillouin zone. In fact, for Lj- ln -C to it is sufficient to focus on the vicinities of the two 
Dirac points, as will be seen below. 

The rest of the calculation is tedious, but straightforward. Let us focus on the 
limit Wi n <C to first. In the vicinity of the K point we write k — K + p and expand: 

Hk = v(p x - ip y ) - a 3 (p x + ip y ) 2 + 0(p 3 ), (10a) 
iF k = F (e x - ie») + F 1 ( Px + i Py )(e x + ie») + 0(p 2 ). (10b) 

The expansion around the K' point can be obtained by flipping e x —* — e x , p x — > 
— p x . Note that going one order beyond the Dirac approximation is necessary, since 
the leading contribution to A4' J ' vanishes. (Indeed, the Dirac Hamiltonian has a 
continuous rotation symmetry, so any third-rank tensor must vanish; the continuous 
rotation symmetry is lifted by the trigonal warping term.) In the nearest-neighbor 
tight-binding model one has 

3 3 2 dt Q 3adt 

v = -t a, a 3 = -t a, F = -3— , F 1 = — tt^t- ( n ) 
2 8 oa 2 oa 

Still, in Eqs. (| 1 a|) . (| 10 &[) we prefer to keep all four parameters v,as, Fq, Fi 

independent, since going beyond the nearest-neighbor tight-binding model would 

invalidate Eq. (fTTj) , while the form of Eqs. (jlQqp , (|106|) is fixed by the symmetrjj§|. 

It is instructive to write down explicitly the expression for the matrix element 

after the integration over the directions of p (we change the integration variable from p 

§ In fact, the Cq v crystal symmetry allows a term ao(p^ m ^ ne expansion of (electron-hole 

asymmetry), which appears already in the second- nearest-neighbor approximation. However, it does 
not contribute to the matrix element calculated here because of its full rotational symmetry. 
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to ek and denote by Ac = L x L y /N — \/27a /A the area per carbon atom): 



\Mn^out V 2MWp h 




Ac CinUJin + Cout^out + Cph<^ph 

Wph u \ 2w ph 2w out 



(12a) 
(126) 



^ r, F0CX3 w out i*l / ^out ^out 1 \ > 
Oout — — ^ — 5 -z V = 1 I , \ LZc ) 

C ph = -2^ + ^-^. (12d) 
tr v 2u in uj ou t 

The integral over ek corresponds to the sum over the intermediate states whose 
energies are 2ek (an electron with the energy — ek in the valence band is promoted to 
the conduction band where its energy is ek)- The poles are bypassed by adding an 
infinitesimal imaginary part ek — ► £k — i0 + , as follows from Eq. (j4|). 

At frequencies 0Ji n not small compared to to the k integral in J\A xxy should be 
evaluated numerically. We focus in the limit w pn <C Wj n . In this limit the overall 
Raman efficiency, i. e., the absolute probability for an incident photon with the 
polarization e m and frequency u>- ln to be scattered in the full solid angle Att with any 
polarization, accompanied by the emission of a 1580 cm -1 optical phonon, is given by 

27rAr I e 2 \ 2 /LU\„a\ 2 



The dimensionless electron-phonon coupling constant is defined as Ar = 
(V27 /Mw p h)(io 1 dto /da) 2 , which coincides with the definition of Ref. [17] in the tight- 
binding model. The dimensionless function f(u)- m /to) is plotted in Fig. [3l In the 
low-energy limit, Lo- m <C to, its value is /(0) = 1, as follows from Eq. ()12aj) . When 
the frequency matches the van Hove singularity at the M point of the first Brillouin 
zone, Win ~ 2to, the intensity diverges, /(w; n /io) VO^inAo — 2) 2 . This divergence 
is cut off at the scale |wi n — 2to\ ~ max{w p h, 7}. However, in the absence of a reliable 
information on electronic relaxation processes in this region of the spectrum, we prefer 
not to study this issue in detail, leaving the divergence as it is. 



3. Discussion 



First, let us see which electron-hole states contribute the most to the G peak intensity. 
In principle, for u>- m -C to one can imagine three situations. They are illustrated Fig. 2] 
where the w electron dispersion is shown schematically, together with the contributing 
states. In the situation (a) these are the states whose energies are close to the half of 
the excitation frequency, the difference ek— Win/2 being small by some parameter. This 
was shown to be the case for the 2D peak, the second-order overtone of the D peak 
at 2700 cm" 1 , which is fully resonant and |ek — u)- m /2\ ~ 7, where 27 <C uj- ln is the 
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in 

Figure 3. The dimensionless function f(uii n /to) appearing in Eq. (1 1 3 D . 



electron inelastic scattering rate jTSJ [T7] . This is also the case for the doubly-resonant 
D peak at 1350 cm -1 , where |ek — uj- m /2\ ~ w p h <C uo- lri [THUSO]- In the situation (b) 
the difference |ek — LO- ln /2\ ~ o^i n itself, and no small parameter enters. Finally, in the 
situation (c) the Raman process is contributed by the whole Brillouin zone, ek ~ to- 

A natural guess for the G peak would be the option (a): indeed, there is a 
contribution from intermediate states that violate energy conservation by the energy 
~ Wph- However, in Ref. |17j it was found that the main contribution to the integral in 
Eq. © with only two terms (|9 q|) . (|9&|) came from the states £k ~ [option (c)]. In the 
present work we have seen that this contribution, in fact, is canceled by other terms, 
missed in Ref. [17], due to the Ward identity. To see that the correct picture for the 
G peak in fact corresponds to Fig. BJc), let us consider doped sample, i. e., with the 
Fermi energy 6f detuned from the Dirac point. Then transitions involving states with 
£k < | e _f" | are simply blocked by the Pauli principle, so the lower limit of the integral 
in Eq. \l2a\ must be set to |ej7-|. Let us analyze the first term in Eq. (|12a[) (the second 
term is almost equal to the first, the third one is smaller by a factor ~ u>^ h /ujf n ): 

oo 

/CinW in C^k _ Cin^in 

\cf\ 

The value of the integral is determined entirely by the ratio |eir|/cJi n , and no other 
energy scales enter, whether small (w p h) or large (to)' That is, the value of the integral 
in Eq. (|12qj) is determined not just by the vicinity of the pole, but by the whole range 
of energies from ek = to ek ~ u)- m . Numerically, when ep is raised from zero to w; n /4 
(half of the distance to the pole), \Ai xxy \ 2 is changed by about 12%. 

In other words, the uncertainty in the energy of the electron-hole states, 
contributing to the process, is of the order of their energy itself (~ 14). By virtue 
of the energy-time uncertainty principle, the duration of the process (the typical 
lifetime of the virtual electron-hole pair) is ~ l/^in- The relevant length scale is thus 
v/u>i n ~ 3 A for v = 10 8 cm/s w 7 eV • A and u)- m = 2 eV. The Raman process giving 
rise to the G peak is thus extremely local in space, involving just a few carbon atoms, 
as it was noted earlier [21] . Its locality can be also understood from the quasiclassical 



2|e F |)+ln 



2M 



2M 



(14) 



Calculation of the Raman G peak intensity in monolayer graphene 

(a) X"-" (b) (c) 




M K 




M K 




M K 



Figure 4. Schematic view of the it electron band structure (solid curves), with 
the states effectively contributing to a Raman process (pink areas)for uj ln <C to: 
(a) states within the narrow region |e k — ajj n /2| <C <^i n ; (b) a broad region of states 
with — Wi n /2| ~ Wi n , but still -C tp; (c) the whole Brillouin zone, ~ to- 



real-space picture of Raman scattering in graphene [20] . This short length scale should 
be contrasted to much longer scales responsible for Raman peaks produced according 
to Fig.il (a), for example, the doubly-resonant D peak at 1350 cm" 1 [2"51 [2^1 W \ 120] . 

Another consequence of the cancellation of the contribution from ek ~ to is the 
frequency dependence of Iq- The standard textbook dependence I oc w 4 (for simplicity 
we assume u- m « w ut) is obtained for systems whose excited states have energies much 
higher than w m [25] ■ The lu 4 dependence, suggested in Ref. [17], was essentially of 
the same origin: it was assumed that the main contribution to the process came from 
states in the bulk of the Brillouin zone, and thus with high energies. As we have seen, 
this conclusion is wrong, and the relevant states have energies ek ~ Win- This remains 
true even at low frequencies, as the electronic spectrum in graphene has no (or almost 
no) gap. As a result, at low frequencies the dependence is oc ojf n , as seen from Eq. (Tl3|) . 
Nevertheless, a recent measurement of Raman peak intensities performed on graphite 
nanocrystallites gives a tof dependence [151 126] . For comparison, we plot Iq from 
Eq. <P~5]) together with two curves oc ojf n with two different coefficients in Fig. [5J 

Let us compare Iq to the intensities of other peaks which have been calculated 
theoretically and measured experimentally. The D peak at 1350 cm -1 is due to 
emission of the scalar A\ phonons with wave vectors near the K point. Due to 
momentum conservation, the D peak is activated by defects or edges. For a defect-free 
graphene flake of the size L a its intensity is given by [20] 

" ' vL a , col + (47) 2 



C\ K I - 



C 2 Lfr.Li 



In- 



(15) 



J D "*»v ( 4 7) 2 

Here uid ~ 0.17 eV and Xk are the frequency and the dimensionless electron-phonon 
coupling constant for the corresponding phonons, v w 7 eV • A is the electron velocity 
(the slope of the Dirac cones), 27 is the electron inelastic scattering rate, L x L y is 
the area of the excitation laser spot, and C is a numerical coefficient which depends 
on the shape of the flake and the character of the edges. We take Ar = 0.03, as 
obtained from the doping dependence of the G peak frequency [6l [7] , and Xk / Ar = 3 
as extracted from the ratio of intensities of the intensities of the 2D peak at 2700 cm -1 
and the 2D' peak at 3250 cm -1 [T7]. Let us assume the inelastic scattering rate to 
be dominated by electron-phonon scattering (which is a lower bound), then it can be 
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Figure 5. The intensity of the G peak as a function of uJi n /to from Eq. d 1 3 1 ) 
for Ap = 0.03, toa/c = 0.0023 (corresponding to to ~ 3.3 eV) (solid curve). 
Dashed and dotted curves correspond to 3 ■ 10 -11 (w in /t ) 4 and 4 ■ 10 -11 (w in /t ) 4 , 
respectively. 



estimated as 7 = (Ap + A/f)wi n /8 [T5]. For a small flake, Eq. (fT3"]) should be weighted 
by the ratio of the area of the flake to the excitation spot area. If we take an ideal 
hexagonal flake with armchair edges (which definitely represents an upper bound for 
Id), then its area is L^a/27/8 and C = 4. For a round flake of the diameter L a with 
atomically rough edges the area is irL^/A and C = 7r 2 /18, i. e., Id is about 10 times 
smaller than for the ideal flake. Taking cj; n = 2eV and using the low- frequency limit of 
Eq. (TTJ), f(cJ- m /t ) = 1, we obtain I D /I G « (300 - 3000 nm)/L„. The experimentally 
found intensity is Id/Ig = (560 nm • eV i ) / (L a ujf n ) [23 [2S], which at u- m = 2 eV gives 
an almost 10 times smaller value than the calculated one for a flake with rough edges. 
We can also consider 2D peak at 2700 cm , whose intensity was calculated to 



For the same values as above we obtain I2d/Ig ~ 150. The experimental value for 
graphene samples on a substrate is typically I2D / Ig ~ 5 0112]. For suspended samples 
the ratio can be I-idJIg ~ 9 [27] and even I2d/Ig ~ 17 [28]. Still, the experimental 
ratio is about 10 times smaller than the calculated one. 

Thus, the theory overestimates both IidIIg an d Id/Ig by about an order of 
magnitude. It seems more reasonable to assume that Eq. (Tl3|) for Iq should be blamed 
for this discrepancy, rather than expressions for Id and Iid- Indeed, as discussed 
above, Id and I^d are determined by just a few material parameters (electron-phonon 
coupling constants, electronic velocity, etc.) which have been calculated by different 
methods and measured in many independent experiments. Moreover, corrections 
to Eq. (p~6j) due to the trigonal warping and electron-hole asymmetry have been 
estimated in Ref. [17] and turned out to be small. At the same time, the parameters 
determining Iq (electronic trigonal warping, electron-phonon coupling anisotropy) are 
not really known, the validity of the tight-binding model at sufficiently high energies 
is questionable, and even the calculated frequency dependence of Iq does not seem to 
agree with the experiment. 



be [BUT] 




(16) 
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Still, what is the main source of such a strong discrepancy? Is it possible that 
Lu- m — 2 eV is already sufficiently close to the singularity at uj- m = 2t , so that the 
dependence of Iq on Lj- ln is noticeably steeper than to 2 , and the value of Iq is higher 
than that predicted by the low-frequency asymptotics? On the one hand, the nearest- 
neighbor tight-binding model with to ~ 3 eV seems to describe the band structure 
reasonably well in the energy range of interest, as is seen from its comparison with the 
results obtained by angle-resolved photoemission spectroscopy [29] . Then the energy 
of electron-hole separation at the M point is 2to « 6eV and u}m/to = 0.6 is quite close 
to the low-energy limit, as seen from Fig. [3] On the other hand, a recent ah initio 
calculation of the band structure of graphite gives the electron-hole separation at the 
M point about 4 eV [30]. Also, the anisotropy of the electron-phonon coupling may 
be stronger than that obtained from the nearest-neighbor tight-binding model |31j . 
In addition, the spectrum of electron-hole excitations can be modified by excitonic 
effects: the Coulomb attraction between the electron and the hole lowers the energy 
of the excitation with respect to its non-interacting value, so that the singularity in 
Id should be at a lower frequency. The author is aware of only one study of excitonic 
effects in the optical response of monolayer graphene; within the Dirac model it was 
shown that excitonic effects result in a weak feature at the border of the electron-hole 
continuum |32j . Their importance for optical excitation of electrons near the M point 
of the first Brillouin zone remains to be studied. Also, a direct measurement of the 
energy corresponding to the singularity by an excitation in the ultraviolet frequency 
range would shed light on this issue. 

4. Conclusions 

We have calculated the frequency-integrated intensity of the Raman G peak in 
monolayer graphene using the nearest-neighbor tight-binding model for electrons. The 
resulting intensity is determined by the trigonal warping of the electronic bands and 
the anisotropy of the electron-phonon coupling, and is proportional to the second 
power of the excitation frequency. Comparison to the intensities of other peaks, which 
have been measured experimentally and calculated theoretically, suggests that the 
present calculation underestimates Iq by about an order of magnitude. 
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